Spectroscopy of Magnetic Excitations in Magnetic Superconductors Using Vortex Motion 
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In magnetic superconductors a moving vortex lattice is accompanied by an ac magnetic field which leads to 
the generation of spin waves. At resonance conditions the dynamics of vortices in magnetic superconductors 
changes drastically, resulting in strong peaks in the dc I-V characteristics at voltages at which the washboard 
frequency of vortex lattice matches the spin wave frequency w s (g), where g are the reciprocal vortex lattice 
vectors. We show that if washboard frequency lies above the magnetic gap, peaks in the I-V characteristics in 
borocarbides and cuprate layered magnetic superconductors are strong enough to be observed over the back- 
ground determined by the quasiparticles. 
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The coexistence of magnetism and superconductivity was 
observed in many crystals, such as RMogSs, RRI14B4, 
RBa 2 Cu307_,5 and (R,A)Cu0 4 _ 5 (A=Sr, Ce) with the tem- 
peratures of magnetic ordering Tm much smaller than the 
superconducting critical temperature T c , and also in borocar- 
bides RT2B2C and ruthenocuprate RuS^GdQ^Os with Tm 
of the same order as T c . Here R is the rare earth element, 
while T=Ni, Ru, Pd, Pt. In such crystals /-electrons of ions 
R give rise to localized magnetic moments, while conduct- 
ing electrons exhibit the Cooper pairing. In all these crystals, 
except HoMogSs and ErRh4B4, magnetic moments order an- 
tiferromagnetically below Tm- Such magnetic ordering co- 
exists with superconductivity without strong interference be- 
cause spin density varies on the scale much smaller than the 
superconducting correlation length and net magnetic moment 
vanishes, for review see Refs.Ml \A X& 



In this Letter we consider interplay between magnetic 
and superconducting excitations in magnetic superconductors, 
particularly interaction between a moving vortex lattice and 
spin waves via the ac magnetic field induced by moving vor- 
tices. The energy transfer from vortices to the magnetic sys- 
tem leads to dissipation which is additional to that caused by 
quasiparticles. This results in strong current peaks in the dc 
I-V characteristics at voltages at which the washboard fre- 
quency of vortex lattice 0] matches the spin wave frequency 
u> s (k) and k matches a reciprocal vortex lattice vector g. 

First we consider slightly anisotropic supercon- 
ductors, i.e. all systems mentioned above except 
SmLai_ x Sr2;Cu04_5 and RuS^GdC^Os crystals, and 
probably also Srri2-. T Ce. Z ;Cu04_,5. The latter are lay- 
ered superconductors with intrinsic Josephson junctions 
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We assume, for simplicity, a uniaxial crystal structure with 
the principal axis along z. The dc magnetic field is applied 
along the z-axis and we assume that the magnetic induction 
B(r), r = x,y, inside the superconductor corresponds to the 
ideal Abrikosov square vortex lattice (such a lattice is real- 
ized in clean borocarbide crystals in field B || c in some field 
intervals Ql). The sublattice magnetization in the case of anti- 
ferromagnetic ordering is assumed to be oriented in the (x, y) 
plane. The dc transport current with the density j is along the 
jy-axis which, due to the Lorenz force, causes motion of the 
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FIG. 1: Vortex lattice moving with the velocity v induces a spa- 
tially periodic ac magnetic field h(x, t) which excites the system of 
magnetic moments shown by red arrows. This additional dissipation 
results in current peaks in the I-V characteristics in magnetic super- 
conductors. 



vortex lattice with the velocity v along the a>axis. 

We use the quasistatic approach assuming that the space 
structure of the magnetic field is the same as in the static vor- 
tex lattice, but the field moves in the same way as vortex lattice 
does. Thus all quantities describing the moving vortex lattice, 
i.e the magnetic field and supercurrents, have the dependence 
on the coordinates and time in the combination (r — vt). In 
the field interval B -C H C 2 the magnetic field should be found 
from the London equations I2L l9l 1 1 Oil 
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curl B = — j 
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47rcurl M, 



J s = 



(1) 

B = curl A, (2) 
curlV0 = ^27n5(r-r„), (3) 



sat V *~*; A 



where j s is the supercurrent, A is the vector potential, </> is 
the phase of the superconducting order parameter, M is the 
local magnetization, $0 is the flux quantum and X± = A^ = 
\ y is the London penetration length for currents in the (x, y) 
plane in the absence of magnetic moments. Further, r n (t) = 
r n(0) + vt are the coordinates of vortices and r„(0) form a 
regular vortex lattice. From Eqs. Q-ifJJi we obtain 



curl curl (B - 47rM) 



A 2 / 



-B 



$0 



S(r - r„). (4) 



To relate the Fourier components of M z (r,t) = M and 
B z (r,t) = B we use the linear response approximation in 
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which supercurrents induce the "external" magnetic field, 

#(k,cj) = S(k,w) -47rM(k,o;), (5) 

acting on the magnetic moments, where M(k, w) = 
x(k, cj)iT(k, w) and x(k, w) = %2z(k, w) is the susceptibil- 
ity of the magnetic system. This approach is valid for the 
magnetization harmonics g x ^ satisfying the condition 



|A/(g, 5 ,«)| 2 /(™) 2 «l. 



(6) 



For an antiferromagnet with two sublattices the magnetic sus- 
ceptibility is given QujllT^I] by 



X(k,w) 



U M UJ s (\t) 



^ 2 (k) 



(7) 



Here ujm = A* 2 «A//(2ft) at fiB <C /cbTm, the density of 
magnetic ions is nu and their magnetic moment is fi, uj s (k) 
is the magnetically active spin wave dispersion renormalized 
by the superconductivity |3J, while v s is the relaxation rate of 
spin waves due to the interaction with phonons. Using Eqs. (0} 
and (0, we obtain for the Fourier components (k = g = 
27r(_Bo/ < l ) o) 1 ^ 2 ("': m i 0); u> = g x v) of the magnetic field 



A 2 , k 2 



1 + 47rx(k, uj) 



B(k, to)=Y, B a S(k - S )5(uj - g x v) 



(8) 

where Bq is is the average induction and n, m are integer. 
From Eq. l|8} we see that magnetic moments renormalize the 
London penetration length so that the effective penetration 
length in magnetic superconductors is given by 01 



A x (k,w) 



A, 



[l + 4 7 r X (k,a;)] 1 /2 



(9) 



Solving Eq. we obtain the Fourier components of the mag- 
netic field B and the "external" field H as 



B(k,w) =B Q 
H(k,u) = B 



8(k- g)6(uj - g x v) 
1 + A 2 ± (g,uj)g 2 ' 
6(k- g)5(uj - g x v) 
l + 47rx(g,w)+A 2 L5 2 



(10) 
(11) 



Thus the moving vortex lattice induces a spatially periodic ac 
"external" magnetic field h(r, t) = H(r, t) — Bq along the z- 
axis characterized by momenta g and washboard frequencies 
oj = vg x . At x = for Aj^ = 1300 A, typical for borocar- 
bides, the amplitude of the main harmonic, n = 1, m = 0, is 
about 20 G. The moving vortex lattice induces also an electric 
field E = [v x B]/c along the current direction. 

When the alternating magnetic field h(r, t) is not parallel 
to the sublattice magnetization, it excites spin waves with mo- 
menta g and frequencies cj s (g) = g • v (this condition of 
resonance holds for any vortex lattice). Assuming that sub- 
lattice magnetizations are almost perpendicular to the applied 



magnetic field, we obtain for the power per unit volume trans- 
mitted from the vortex lattice to the magnetic system the ex- 
pression fl2ll 



V 



M 



M(r,t) 



dh(r,t) 

dt 



2g x v\h( S , 5^)| 2 Im[ X (g, g x v)}, (12) 

g 

where angular brackets denote time and space average. 

Now we are in a position to find the velocity of the vor- 
tex lattice at a given transport current density j. For that we 
equate the power per unit volume performed by the battery, 
jE, to the sum of the power dissipated by quasiparticles, rjv 2 , 
and that transmitted to the magnetic system, Vm- Here r\ is 
the viscous drag coefficient due to quasiparticles in normal 
vortex cores. It is given by the Bardeen-Stephen expression 
r/ = BoH* 2 o- n /c 2 , where a n is the normal state conductivity, 
H* 2 = $o/(27r£j_) is the orbital upper critical field and £j_ 
is the superconducting correlation length in the direction per- 
pendicular to the applied magnetic field. Taking into account 
that E = vBq/c and u = vg x = cEg x /Bo, we find v and fi- 
nally j-E (i.e., I-V) characteristics in the intervals of E, where 
inequality Eq. (jfji is fulfilled: 

• /pr\ = ° 2 V F V" 2g x cB a lm[x(g,cEg x /B )} 

B ° ^ + ^x{^cEg x /B )+^\\\^ 

(13) 

From this equation we see that the current density as a func- 
tion of E has peaks corresponding to resonances between 
the ac magnetic field and spin waves, i.e when u)(n,m) = 
2irv{B /<S> ) 1/2 n. 

Let us discuss the behavior of j(E) near resonances. We 
introduce the frequency deviation Aui = uj s (g) — oj such 
that v s <c Aw <C u> s (n,m). Then we obtain x(g,w) w 
lom/(2Auj) and Im[x(g, uj)] « ljmv s /(2Alu) 2 . We con- 
sider the interval of frequency deviations Aui where X]_g 2 S> 
47rx(g, uj). In this interval we estimate 



M(g,u>) ^ 

/UHm 167r 2 (n 2 + m 2 )\ 2 L h/S.uj 



(14) 



Due to the condition Eq. l|6) our approach is valid for HAui > 
/i<I>o/(47rA^) 2 . The ratio of the additional current caused by 
spin waves over the current background is given as 



Aj(n,m) ujmKs^qBq 



8tt 2 w7?(Aw) 2 A| (n 2 + to 2 



(15) 



In the frequency interval HAuj > /i$o/(47rA^) 2 we obtain 
the inequality 



Aj(n,m) 16ir 2 hriMV s Bo n 2 



uji]^o {n 2 + ni 2 ) 



2\2 ' 



(16) 



In magnetic insulators v s is typically of order 10 6 s 1 . One 
can anticipate the same value in magnetic superconducting 
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crystals, as conducting electrons are gapped. For H0M2B2C, 
taking H* 2 w 10 T, n M = 10 22 cm" 3 , a n = 10 5 
(ohm- cm) -1 at u> = 10 10 s _1 we derive Aj(n,m)/j < 
0.8n 2 /(n 2 + m 2 ) 2 . Thus, the peak n = l,m = is observ- 
able even in the frequency interval where our linear response 
approach is valid. Here the magnetic system deviates only 
slightly from equilibrium as energy is transformed further to 
phonon bath. 

Closer to the resonance the linear response approach breaks 
down. Here the dominant contribution to dissipation comes 
from generation of spin waves by vortices which leads to a 
strong deviation of the magnetic system from equilibrium. For 
quantitative description of the j-E characteristics close to res- 
onances the full dynamic approach for vortices and magnetic 
system is necessary. 

Based on Eq. (II 3i we see that measurements of the I-V 
characteristics at different magnetic fields and currents may 
provide information on the spin wave dispersion w s (g). The 
washboard frequency oj and the reciprocal vortex lattice vec- 
tors g may be changed independently by varying Bq and j, 
but an important question is what are limitations on the vari- 
ations of the magnetic field and the current density. Mo- 
mentum k ~ 2-7r(i?o/$o) 1 / 2 is of order 10 6 cm -1 in fields 
Bq < 1 T and increases as one approaches H c2 , but then 
harmonic amplitudes h(g, oj) drop. Limitations on frequency 
are due to limitations on the current density, which should be 
lower than the depairing current density, and also should not 
lead to excessive heating. From Eq. (1131 . to reach frequency 
uj one needs current density j(ui) > a n uiH* 2 /cg x and the 
electric field E(ui) = ujBo/cg x . For Bq = 1 T we obtain 
j(cu) as 10 7 n _1 (7iu;/lK) A/cm 2 if the lowest harmonics are 
used, while for higher harmonics higher frequencies may be 
reached. The depairing current density for borocarbides is of 
order 10 7 A/cm 2 and thus spin waves with energies Huj < 1 K 
may be probed without strong suppression of superconducting 
properties by the transport current. For the dissipation power 
per unit volume, Vdis = jE > a n uj 2 ^oH* 2 / '47r 2 c 2 , we esti- 
mate Vdis ~ 10 8 n~ 2 (frcj/lK) 2 W/cm 3 . To diminish heating 
the pulse technique may be used, as in I-V measurements by 
KunchurQ. 

As only the low energy part of the spin wave spectrum may 
be probed by I-V measurements, the important question is 
what is the minimum energy of spin waves. For a Neel anti- 
ferromagnetic state there is always magnetic gap in the spec- 
trum due to the magnetic anisotropy slightly renormalized by 
superconducting screening of RKKY and dipole-dipole inter- 
actions 0]. Neither experimental nor theoretical information 
on the strength of magnetic anisotropy or the structure of ex- 
citations in these materials is available so far. Thus we cannot 
predict yet whether resonance conditions for the lowest har- 
monics will be fulfilled in borocarbides. However, we can 
anticipate that higher harmonics will be effective in the case 
of weak pinning. 

Next we discuss highly anisotropic layered crystals like 
SmLai_ x Sr 2; Cu04_ ( 5. This material is especially interesting 
because its T* structure leads to a two-dimensional character 



of the magnetic system. Here magnetic Sm2C»2 and nonmag- 
netic La2-. T Sr T 02-a layers alternate in the barriers between 
the superconducting CuC>2 layers. The Josephson nature of 
the interlayer coupling in this crystal has been confirmed by 
observation of the double Josephson plasma resonance stem- 
ming from two layers in a unit cell J^. The specific heat mea- 
surements |0] show that magnetic ordering is absent down to 
a temperature of 0.3 K and a magnetic gap, if any, lies below 
0.3 K. They reveal also a broad peak near the temperature 1 K 
and the height of this peak indicates the presence of competing 
interactions that might be described by the two-dimensional 
J1-J2 Heisenberg model with J2/J1 > 0.4 HIS]. Such 
a model has very complex dynamics and contains a variety 
of transitions down to zero temperature, making it an ideal 
testing ground for the theory of quantum phase transitions. 
The most interesting part of the phase diagram is in the re- 
gion 0.4 < J2/J1 < 0.55, where a gapped phase without 
magnetic ordering is likely to be taking place. However, its 
characterization has been one of the most intriguing puzzles 
of the physics of strongly correlated systems 1 16]. 

If the magnetic field is applied perpendicular to the layers 
(along the c-axis), it induces pancake vortices which do not 
form a regular lattice in magnetic fields above 20 G as they 
order along the c-axis only due to weak Josephson and mag- 
netic interactions fTvIl . This makes excitation of spin waves 
ineffective by moving vortex lattice induced by a perpendic- 
ular magnetic field. When a magnetic field is applied paral- 
lel to the layers (in the a6-plane, along the y-axis), the situ- 
ation is drastically different, because now Josephson vortices 
H d |H I2H |2l are induced. In high fields they form a 
lattice which is quite regular in the x-direction (parallel to the 
layers). Josephson vortices do not have normal cores and so 
only thermally induced quasiparticles (or those near the nodes 
in the case of d-wave pairing) cause dissipation. A weak in- 
terlayer tunneling transport current, which leads to vortex mo- 
tion in the x-direction, cannot destroy superconductivity and 
produces much less heating than in the case of isotropic or 
weakly anisotropic superconductors. 

The distribution of the magnetic field B(r) inside intrinsic 
Josephson junctions is described by coupled finite-difference 
differential equations for the phase difference ip n and for the 
magnetic field B n inside the junction n between layers n and 
n + 1 I20II2I. Accounting for the magnetization M n of ions 
inside intrinsic Josephson junction n we obtain equations for 
the dimensionless variables cp n , b n = B n 2ir\ a t,\ c /Q , m n = 
M n 2ir\ ab \ c / $ and h n = b n - 4irm n : 



d 2 ip n 
dr 2 



diPn 

c dr 

bn dtp n 

I 2 du 



sin ip r , 



= 0, 



+ Vab 



_ dhn 
du 
d_ ( df n 
dr \ du 



(17) 



0. 



where u = x/Xj, r = tuj p , Xj = js, s is the interlayer 
distance, 7 = X c /X a b is the anisotropy ratio, A c and X a b 
are the London penetration lengths for currents along the c- 
axis and in the a6-plane, respectively, uj p = cj (Acy^) is the 
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Josephson frequency, e c is the dielectric function along the 
c-axis, v c = 4ttct c / (u> p e c ), v ah = Ai:a ab /(^ 2 e c u; p ), a c and 
a a b are quasiparticle conductivities along the c-axis and in the 
a6-plane, respectively. Using linear response approximation, 
m„ = X^n/(1 + 47rx), where \ = Xyy> we see that h n 
h, : I -+- 4ttx) satisfies the same equations as b n at \ = 0, 
but with the renormalized parameter £~ 2 = (1 + 4irx)£~ 2 - 
For SmLai„ 2 ;Sr :1 ;Cu04_i we estimate Cr 2 <C 1 because 
lu m ps 1.8- 10 s s- 1 , € 2 ps 2-10 4 at^ = 0.8/i B , n M = 5-10 21 
cm" 3 and A a6 ps 2000 A. 

In the following we consider large enough fields B > Bj = 
$o/(27rsAj). Then the Josephson vortices fill all intrinsic 
junctions, overlap strongly and form a regular trian gula r lat- 
tice GHHEHIIl. (An illustration is given in Ref.[ll) For 
SmLai_ x Sr 2; Cu04_a we have 7 ps 500, oj p ps 10 12 s _1 and 
£?,/ ps 0.5 T. In a Josephson system the washboard frequency 
is the Josephson frequency co = uj.j = 2eV/H, where V is the 
voltage between neighboring layers. For a triangular lattice at 
frequencies and the magnetic fields satisfying the conditions 
I 2 S> (1 + 4-7tx) and \2u> — 6| > 1, where Co = uj/uj p and 
b = Bq/ Bj, the solution of Eqs. dl7> has the form 



tp n (u, t) ps Cut — bu + nn + 



4 sin(a)r — bu + 7m) 
4w 2 - b 2 : 



h n (u, t) ps — /io cos(wr — bit + 7m), ft-c 



4tD 2 - fe 2 ' 



where we neglected 2^ and !/„&. We estimate h = 
/i $o/(27rA 2 b 7) « 0.16 G at lu = 0.1u p and B = Bj. Near 
the Eck resonance, 2Co ps 6, the amplitude of the magnetic 
field h is larger. For the reciprocal lattice vector we have 
g = (27rs J B/$ ,0,7r/s). So g x = 1/Xj ps 10 4 cm" 1 at 

B = Bj. 

Assuming that sublattice magnetization is almost perpen- 
dicular to the applied magnetic field or that magnetic ordering 
is absent we obtain for the I-V characteristics 



V esh 2 



C e ff- 



-Im 



Xyy I S 



2eV 



(18) 



where a e s = a c + 2a a bB 2 /('yB) 2 describes dissipation 
due to quasiparticles. At resonance uj = w s (g), we es- 
timate Aj/j ps 2n 2 c 2 s 2 h 2 uJM/{uj<J c ffVs®o)- Estimating 
a c ~ 10- 3 (ficm)- 1 and a ab ~ 4- 10 4 (ficm)" 1 as in BSCCO 
and taking s ~ 12 A as in Sm2-xCe x Cu04-,5, we obtain 
Aj/j ps 4 and |M(g, wj)|/(/in A /) ps 0.3 at uj = 10 12 s" 1 
and i? = Bj and bigger values near the Eck resonance. Cer- 
tainly, such frequencies are sufficient to probe almost com- 
plete spectrum in SmLai_ x Sr 2 ,Cu04_5. 

In conclusion, we propose to probe low-frequency mag- 
netic excitations in magnetic superconductors by measuring 
I-V characteristics in the mixed state with a moving vortex 
lattice. Coupling of such a lattice to magnetic moments is due 
to an ac magnetic field which is inherent to vortex motion. The 
energy interval of spin waves which can be probed in isotropic 



and moderately anisotropic superconductors is limited by the 
depairing current and heating. If spin wave energies fall in 
this interval, they affect the vortex motion strongly and should 
be easily seen in the I-V characteristics as current peaks at 
corresponding voltages. Such an effect may be observed in 
borocarbides if they have spin waves with energies below 1 
K. For highly anisotropic layered superconductors in parallel 
magnetic fields, higher spin wave energies may be probed by 
use of moving Josephson vortices. This is sufficient to study 
almost complete spin wave spectrum in SmLai-zSr^CuC^-i 
with exotic magnetic ordering. 
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